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One of the most convenient methods to obtain information about the energy distribution function 
of electrons in conducting materials is the measurement of the energy resolved current j{ui) in field 
emission (FE) experiments. Its high energy tail j>(u>) (above the Fermi edge) contains invaluable 
information about the nature of the electron-electron interactions inside the emitter. Thus far, j > (co) 
has been calculated to second order in the tunnelling probability, and it turns out to be divergent 
toward the Fermi edge for a wide variety of emitters. The extraction of the correlation properties 
from real experiments can potentially be obscured by the eventually more divergent contributions of 
higher orders as well as by thermal smearing around Ep. We present an analysis of both factors and 
make predictions for the energy window where only the second order tunnelling events dominate the 
behaviour of j > (u}). We apply our results to the FE from Luttinger liquids and single- wall carbon 
nanotubes. 
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I. INTRODUCTION 

According to the original idea of Fowler and 
Nordheimi, the cold emission of electrons from metal- 
lic electrodes can be considered as tunnelling of particles 
from the conduction band into the vacuum through a tri- 
angular barrier. The precise form of the latter is deter- 
mined by the work function Wa, by the applied electric 
field F as well as by the geometry of the emitter tip, 
and is characterised by the effective transmission coeffi- 
cient D(u>) for the electrons with energy ojS. The most 
important measurable quantities are the energy resolved 
current j(u>) and the total current J = J dioj{uj). The 
currents drawn from the emitters are usually very small 
so that j(uj) is proportional to the equilibrium energy 
distribution function of electrons in the vicinity of the 
tip, np{u>), and to 



j(uj) ~ D{uj)n F (uj) 



(1) 



Therefore, at T = no electrons are allowed with ener- 
gies above the Fermi edge Ep. However, in 1970, Lea 
and Gomer succeeded in measuring j(ui) to a very high 
degree of accuracy. They found significant contributions 
above Ep£. This effect can be attributed to electron- 
electron interactions inside the emitter. The relevant 
process is usually referred to as "secondary tunnelling" 2 , 
during which an electron with energy below the Fermi 
edge tunnels out of the emitter leaving behind a hot hole. 
In the second stage, this hole scatters inelastically, pro- 
ducing a secondary electron which can obtain an energy 
above Ep. Due to its higher energy, this electron has an 
increased tunnelling probability and will contribute to a 
"secondary current" j>(u)) above Ep. Since j>(oS) comes 
about as a result of interactions, its measurement reveals 
valuable information about the nature of correlations in- 
side the emitter. 

Two different types of rigorous calculations of j > (u) 



have been attempted so far. The approach o£a used ki- 
netic equations. Ivm the high-energy tails in j(to) were 
calculated using non-equilibrium perturbation theory in 
the tunnelling probability. One of the principal predic- 
tions is the asymptotic behaviour toward the Fermi edge, 
which turns out to be divergent. Both approaches take 
into account only tunnelling processes of lowest order. 
To the best of our knowledge thus far no investigations 
of processes of higher orders have been undertaken de- 
spite the fact that they can potentially be even more di- 
vergent around Ep and therefore obscure the predicted 
effects. With this work we would like to address this 
problem and analyse the influence of higher order tun- 
nelling processes on j> (w) . 

Thermal broadening of the Fermi edge can also pre- 
vail over the contributions of lowest order, making their 
measurement difficult if not impossible. While in con- 
ventional Fermi liquid systems this broadening is rather 
simple to estimate, in many interacting systems it has not 
yet been analysed. One of the most important of these 
is the universality class of ID correlated electrons or the 
so-called Luttinger liquids (LLs), which are e.g. realised 
in single- wall carbon nanotubes (SWNTs), as has been 
shown both theoretically and experimentally 6 ' 7 ' 8,9 . 

Recently, the advancement in manufacturing enabled 
the construction of the first flat panel display which em- 
ploys carbon nanotubes as emittersAfiiAi. Refinement and 
further development of this technology requires detailed 
knowledge of their emission spectra. This is the reason 
why it is important to fully understand the field emission 
properties of SWNTs. Thus far, theoretical activities in 
this field have been rather moderate. In this work, we 
would like to analyse the effects of tunnelling processes of 
higher orders from LLs along with the effects of thermal 
broadening of the Fermi edge on the high-energy tail of 
j> (w) and apply the results to FE from SWNTs. 
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FIG. 1: Schematic view of the system under consideration 



II. NEXT-TO-LEADING ORDER IN 
TUNNELLING: DIAGRAMMATIC APPROACH 

We model the field emission set-up by a heterogeneous 
tunnelling junction (Fig.^l. The left contact, which will 
be referred to as "host" , represents the emitter. The 
chemical potential of this system is set to fj.^ = 0. We 
shall assume that tunnelling only takes place precisely at 
the interface between the conductors. The right contact 
will be referred to as "lead" , and is represented by a non- 
interacting Fermi gas with chemical potential = —eV. 
As during field emission electrons tunnel into vacuum, 
this potential plays the role of a cut-off and has to be 
set to the conduction band width or the characteristic 
energy scale entering D(uS), whichever is larger. The 
transmission coefficient D{uj) of the triangular barrier en- 
countered in field emission depends exponentially on the 
energy and the work function W a- We restrict ourselves 
to the low temperature limit, where interesting processes 
will only involve electrons quite close to the Fermi edge. 
That is why we replace this tunnelling probability by a 
constant, 7 2 ~ -D(O), where 7 is the energy independent 
tunnelling amplitude (we set Ep = throughout). The 
model Hamiltonian is then given by 

H = H[i/>] + H[c] + 7 ['0 t (O)c(O) + V(0)c f (0)], (2) 

consisting of three parts. H[tfj] and H[c] are the Hamil- 
tonians of the host and the lead, respectively, and the 
third term describes the tunnelling events. 

Due to the finite applied voltage, the system under 
consideration is in a non-equilibrium state. The most 
convenient method of describing such systems is the 
Keldysh formalism^, where we shall be taking iij nt = 
7[?/' t (0)c(0)+V(0)c t (0)] as the perturbation of the Hamil- 
tonian Hq = H[ip] + H[c\. We shall use the notation 
employed ini^, where the Green's functions are defined 
by 

g—(t,t') = -i(T[i;(t)^(t')Sc])o, 
g-+(t,t') = i(^(t'W(t)S c )o, 
g + -(t,t') = -i{m^(t')Sc)o, 
g ++ (t,t>) = -i(f[^(t^(t')S c }) . (3) 



T and T are time and anti-time ordering operators, re- 
spectively, and (• • -)o stands for the averaging with re- 
spect to the ground state of H . We use g(t, t') to denote 
Green's functions of particles in the host, whereas G(t, t 1 ) 
refers to Green's functions in the lead. The operator Sc 
is defined by the perturbation, 

S c = T c cxp j- l -l j dt [tf (t)c(t) + (*)] J • (4) 

As shown inil the energy resolved current for a certain 
energy u > is proportional to the local Green's func- 
tion g *~(uj) in the host. Our primary task is to calculate 
the function g ^(o-i). For this purpose we shall perform 
a perturbative expansion of Sc in the tunnelling ampli- 
tude 7, so that g~ + (u) = g ( ^ + (aj)+g 2 ~ + (cj)+ i gj + (a;) + . . . 
where g^ ~ 7 2 and 34 ~ 7 4 . As the Fermi edge is still 
sharp even in interacting systems, the zeroth order does 
not contribute above Ep. According to&, the second or- 
der correction to the Green's function can be brought to 
the form, 

2 p f poo 

fc+M = ^jdte^jj dhdt 2 (5) 
x ^ (ij)K ij (t,t 1 ,t 2 )G ij (t 1 ,t 2 ), 

ij=± 

where ify are four-point correlators of V'-operators which 
correspond to rectangular diagrams with insertions^*. It 
has been shown that only 

K + .(t,h,h) = -(f{^ + (0)V t (ti)}T{V'(t)V'(*2)})(6) 

yields a non-zero contribution above w = 0, which di- 
verges toward the Fermi edge. 

We can use the same approach that led to |[5J| and 
expand to the next-to-leading order in 7. It turns out to 
be given by 

91 + M = |£ / dt ^ t/h ////_" dtl ' ' ' du (7) 

x K ^ kl (*' *i. *a. *3, U)G ij (h , t 2 )G M (t 3 , U) , 

ijkl 

which schematically corresponds to a sum over 16 differ- 
ent hexagon-shaped diagrams, s. Fig. [21 each containing 
two insertions of the lead Green's functions and a six- 
point correlator which is similar to the four-point 
correlators in the second-order case. The structure of 
the diagrams is not surprising as in this order two suc- 
cessive tunnelling processes are allowed. Since we are 
only interested in the high-energy contributions above 
Ep, we first have to identify the diagrams which may 
contribute in this situation. As in the case of lowest or- 
der contributions, the best way to identify the relevant 
constellations without explicit calculation is to take ad- 
vantage of the spectral representation. To that end we 
have to be able to perform the t- integration in Q. This 
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FIG. 2: Fourth order diagram: The hexagon corresponds 
to a six-point correlation function of ip operators and the in- 
dices i,j,k,l = ± denote the respective part of the Keldysh 
contour. Solid lines depict Green's functions in the host, 
dashed lines stand for Green's function is the lead wire. 
Crosses stand for the tunnelling vertices. 

can be accomplished by insertion of an appropriate num- 
ber of complete sets of eigenstates, thereby extracting the 
t-dcpendcncc. 

In order to illustrate the procedure, we shall deter- 
mine the contribution of the K correlation function 

where the time indices t\ to £4 are taken on the time- 
ordered "— " shoulder of the Keldysh contour. Hence, we 
wish to calculate 

<7 4 -+H = ^ J dte^' h JJJJdti---dU (8) 

The time development of the ip(t) and c(t) operators in 
the interaction picture is governed by the unperturbed 
Hamiltonians H[ip] and H[c], respectively. Although due 
to the interactions in the host, we cannot diagonalisc 
H[4>], we still know that there is a complete set of eigen- 
states satisfying H [V>]|^) = E v \v) and ^„ \ v){v\ = 1 with 
positive eigenenergies E v . Inserting a complete set into 
© yields 

x {Q\e- iH ^ h ^e iH ^ h \v)A{v), 
where A(y) is independent of t, 

A{v) = JJJJ dtl '" dti ( 10 ) 
x (u\T[^(0)4^ 2 4^\0)(T[cxclc 3 ct]} . 
Hence, the t integration can be performed and leads to 

, 94 -+H = i^E/ dte^e iE ^(^\v)A{v) 
= %H^ + E v m^\u)A{u). 

As all energies E v are positive, we conclude that this 
function does not yield a contribution above to = 0. Sim- 
ilar calculations can be performed for all other correla- 
tion functions. We find out, that e.g. in second order, the 



only contributing diagram is described by the correlation 
function, 

k+- = -(f[^(o)4m(t)^})> (11) 

which contributes for to < eV, the upper threshold for 
second order processes. On the other hand, in fourth 
order there are only three relevant contributions which 
are given by 

K+ = (f[^(0)4]Tm)Mi^}) (r (12) 

K+-+- = (f[4(o)ipi4m(t)i>2iP4) o , (13) 

K+++- = (f[4(0)4^l}Tmt)^]) Q - (14) 

Two of these correlation functions (K^ and K ++ ^ ) 

contribute up to w = eV and have to be counted twice 
due to reasons of symmetry. The correlation function 

1 contributes for to < 2eV which is the upper 

threshold for fourth order processes. In all of these con- 
tributions at least one of the insertions is of the type 
encountered in the lowest order. We know that in the 
lowest order (the rectangular diagram o£i4) the purely 
time or anti-time ordered insertions do not lead to con- 
tributions above to = 0. Therefore, in the terms where 
only one of the insertions is of the type G + ~(uj), one 
can regard the corresponding diagram as an effectively 
rectangular one with respect to high-energy tails. 



III. FIELD EMISSION FROM LUTTINGER 
LIQUIDS AND SWNTS. 

The usual way to proceed after the relevant diagrams 
have been identified is to perform a perturbative expan- 
sion of all contributions in the interaction constant. How- 
ever, a perturbative breakdown of the six-point correla- 
tion functions (hexagons of Fig. [2J in non-equilibrium 
involves a very large set of diagrams. Performing the 
same calculation for a LL system does not have that 
disadvantage. The reason for that is the fact that in 
the bosonization representation of an LL any possible 
fermionic correlation function can be computed for any 
interaction strength without resorting to approximation 
schemesi^iSiii. The perturbative results are then ex- 
pected to be recovered after taking the limit of vanishing 
interactions. On the other hand, as we are showing later, 
the results for an LL only have to be slightly modified in 
order to be applied to the experimentally relevant FE 
from SWNTs. 

Our main goal is the calculation of the kernels Jf^H . 
We assume the emitter to be infinite extending from x = 
—00 to x — 0. Hence, we are dealing with a semi-infinite 
wire with an open boundary condition (OBC) at its end. 
Using OBC bosonizationi^, exactly at the emitter tip 
(x — 0), where all tunnelling events take place, we have 
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for the electron field operator 
ip{x = 0,t) 



'2na 



(15) 



where a is an infinitesimal cut-off length which can be in- 
terpreted as the lattice constant of the underlying lattice 
model. We omit the "ladder operator" (particle number 
changing operator) as we are later going to the limit of 
infinitely long wire L — > ooiiiS. The bosonic phase field 
4> entering (|15|l is given by 



(x = 0,f) 



1 x / 7T 

2^^ V Lq 

9>0 V 1 



e ^-a n + e l " vt at)e~ aq/2 . 



-iqvt 



(16) 

In the above equations, g is a dimensionless quantity 
measuring the interaction strength^, 



Up 
irhvF 



-1/2 



(17) 



where Uq is the Fourier transform of the screened 
electron-electron interaction potential U(x — y) at q = 
and vf is the Fermi velocity, v = vf/q- 

By using the Bakcr-Hausdorff formula, we can calcu- 
late arbitrary correlation functions of the type 



= <«!••• a N ), 



(18) 



where N is an even number and oti are fermion creation 
or annihilation operators, ^{ti) or ip(ti). We arrive at 
the result, 



2y 



F(0) N / : 



N-l N 



F{t m - t n ) 



m—1 n— m+1 



i/s 



where 



F(t) = 1 



-ieo(t— iS) 



(,19) 



(20) 



The parameters eo and 5 are defined by eo = K / Lv and 
8 = a/v. Furthermore Pi = +1 [j3i = —1] for oti = tpiU) 
[on = ^ (*.)]. 

Using l|19|) and taking into account the time orderings, 
we obtain for the rectangular diagram^ 



9 2 + M 



7 



dte luJt/R 



dtidt 2 



(2irah) 2 

x sgn(-ti)sgn(i - t 2 )G + ~ (ti,t 2 ) 

■ F(0)*F(-\h\)F(\t-t 2 \) 
F{-t)F(-t 2 )F(t 1 -t)F(t 1 -t 2 ) 



i/s 



(21) 



Assuming the lead to be a Fermi gas with chemical po- 
tential fiR ; = — eV, the corresponding Green's function 
becomes 



G + -{t u t 2 ) 



Kv e ieV(t 1 -t 2 )/h 

2ir ti — ti — iS ' 



(22) 



with an infinitesimal parameter S. For reasons of sim- 
plicity, we have assumed a constant density of states 
v in the lead. Furthermore, we approximate F(t) by 
F(t) w ieo(t — iS), which is justified for L — * oo. These 
approximations allow the calculation of g^~ + (cj), which 
has been performed ini^. However, the approach oi^ 
cannot be adapted to the fourth order. Nevertheless, 
as 7 is very small, g7 + (w) can become important only 
in the immediate vicinity of the Fermi energy and even 
that only when it turns out to be more divergent than the 
lowest order term. Therefore, it is sufficient to find the 
lu — > asymptotics and we require eV/ui — » oo. In that 
case, the Green's function G" 1 can be shown to become 
a ^-function, 



G + -(t 1 ,t 2 ) -► -ifw5{tx-t 2 ). 



(23) 



Using this representation, the high energy particle num- 
ber n 2 {Lo) = —ig2 + (u))/H can readily be calculated and 
one immediately obtains 



n 2 (uj) 
where lu c = 7/ a and 



1/9-2 



(2 



— (-) 

7t) 2 \vh) 



1/9 



(24) 



(25) 



k 2 contains a numerical prcfactor which is given by the 
dimensionless convergent integral 



k 2 = H) 1/9 lim 

<5^0 



dye ly / d?/isgri(-j/i)sgn(j/ - 



{-\yi\){\vvi\) 



{-y - iS){~Vi - iS)(yi - y-iS) 



1/9 



(26) 



Within the approximation eV/to — > 00, the calculation of 
the fourth order contributions can be performed likewise. 
We start with the determination of the contribution of 

the if-i 1 correlator. In analogy with J2J, we obtain 

a similar expression but with two insertions. According 
to (|23|l we obtain two S functions which allow to perform 
two of the four integrations. Switching to dimensionless 
quantities, we obtain 



.94 + M 



ihv 1/7 . / 7 
-K A — t 



1/S 



_7_ 

aw 



3-1/9 



(2tt) 3 2 '" 4 \vh 
where the numerical factor k' A is given by 



(27) 



K = {-i) 1,g J dye* J J d Vl dy 3 (28) 
sgn(-yi)sgn(-j/ 3 )sgn(y - 2/i)sgn(y - y 3 ) 



(-y-iSy/9 
1 12/3 1 - 2/3 1 



1/3 



(-2/1 - i5)(-y 3 - i<5)(t/i - y - id) 

[j/ — Z/x I la/ — 2/3 1 -2/3 1 

(2/1 - 2/3 - iS)(y 3 - 2/ - i5)(#3 - 2/1 ~ *<$) 



1/9 
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The other diagrams, (ijkl) = (H ) and (ijkl) = (++ 

H — ), can be treated in the same way with the exception 
that the time ordering combinatorics is significantly more 
complicated. However, we found that all three diagrams 
yield the same leading power-law energy dependence, 



714(0;) 



1/9-3 



with 



v ^7 / 7 \ 1/9 



(29) 



(30) 



For the (ijkl) = (H 1 — ) case the corresponding k' 4 is 

given by (|28|) . Numerical evaluation of k\ yields a value 
around unity. Due to the different time ordering com- 
binatorics of the six-point correlation functions, the nu- 
merical prefactors k'l and k'l' for the other two diagrams 
cannot be written in a compact form similar to (|28|l . As 
the other two diagrams are related to those of the second 
order in tunnelling (see Section [HJ, k'l and k'l' have to 
be of the same order of magnitude as k'^. 

Equation (|29|l is the central result of this paper. The 
contribution of next-to-leading order processes turns out 
to be even more divergent than the leading order in case 
of weak interactions, g > 1/2. However, we expect the 
overall prefactor to be negative as in the fourth order 
in tunnelling the recombination of the secondary holes 
and secondary electrons, stemming from different pri- 
mary processes, is allowed, effectively suppressing the 
high-energy tails in j > (u>). Taking the limit g — > 1 we 
expect to obtain the asymptotics of the next-to-lcading 
term for a Fermi liquid, m(uS) ~ 1/uj 2 . Similarly to the 
leading order tunnelling results, at g = 1 j>(u>) is identi- 
cally zero due to a vanishing prefactor k^. That can be 
shown analytically. 

So far, all calculations have been made for the case of a 
spinless Luttinger liquid. However, carbon nanotubes are 
known to be described in terms of a four-channel LL&L2S. 
The four channels are responsible for the charge-flavour 
(<f) c -), total spin (</> s -), spin-flavour (4> s +) and total 
charge (<t> c +) excitations. The three former arc non- 
interacting (g = 1) whereas the latter possesses the Lut- 
tinger parameter g c+ = (1 + iVo/nhvp) -1 ^ 2 , where Vo 
is the k — component of the Fourier transform of the 
screened interaction potential. It can further be shown 
that the field operator ip can be written in terms of Bosc 
fields as 



2na 



exp \ji 



/v^T + c - + S + + 4>b-)] ■ (3i) 



As different species of <f> operators commute, all corre- 
lation functions factorise. Therefore, the results for the 
one-channel LL can be applied to SWNTs by the simple 
replacement 



9cl 



(32) 



An estimation of the typical interaction parameter g can 
be done using the following formula, which relates the ge- 
ometry of the nanotube to the interaction constant g c +^, 



9c+ = 1 



irnhvF 



[ln(L/27ri?) + 0.51] 



-1/2 



(33) 



where k is the dielectric constant, L is the length of the 
nanotube and R its radius. Assuming typical values we 
arrive at an interaction constant g c+ m 0.2 which leads 
to an effective g m 0.5. 



IV. VISIBILITY OF HIGH ENERGY TAILS 

Thus far, we have discussed the zero temperature case, 
where a non-trivial high-energy tail in the energy resolved 
current is produced by the interaction effects. However, 
such a tail can trivially emerge as a result of the ther- 
mal broadening of the Fermi edge. In this Section, we 
shall address the problem of distinguishing the correla- 
tion effects from those of finite temperature restricting 
ourselves to the model of spinless LL. 

We first calculate the width of the broadening as a 
function of interaction constant and temperature. Obvi- 
ously, the main contribution to j>(u>) simply stems from 
the zeroth order particle density for a Luttinger liquid, 
which can be brought to the formSi 



ht(uj) 



1 

ha 



dte iujt ' n - 



Tig 
hilv 



sin 



-^(-ivt + a) 



(34) 



where X = 1/g. In order to perform this integral, is is 
convenient to express the sine function in the denomina- 
tor as an exponential and to rewrite the exponent A in 
the denominator with the help of the binomial formula. 
Then, the integral can be performed and we obtain a 
representation of n(uS) in terms of the hypcrgcomctric 
function F— , 



tit(o;) 



sin A (^f )(2i) A 

h0V JaM + (-1)*JJM (35) 



ha 



with 



hi") 



hj3 



7rA — i/3u> 



(36) 



x m *_!^,i + *_!^, e -^/c**o 

'2 2?r ' 2 2tt ' 



This function only depends on the combination (3uj rather 
than on ui alone. Consequently, we can deduce that 
the width of the thermal broadening will be of the form 
o>c = ksTf(\) with a function /(A) independent of the 
temperature. 

Using an appropriate expression for the hypergeomet- 
ric function in the case A = 1 and using the expan- 
sion for a — > 0, the Fermi distribution function of the 
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non-interacting case can easily be recovered. For the in- 
teracting case, i.e. for A > 1, we need the following 
expansion™, 



F(x,y,z,e^) 



(37) 



T{z)T{x + y-z) 

T(x)T(y) 
( 1 _ e i< P y-x-v [i + 0( e ^ _ 1)] 

T(z - x)T(z - y) 



Neglecting terms of order 0(e l(p — 1) which is appropriate 
in the limit of small a, after some algebraic transforma- 
tions we arrive at 

^ = T^wS V-AW.A) (38) 

p -/3u>/2 



[cosh(/3u;) - cos(ttA)] T(l 



2 2ir • 



Having established this relation for finite tempera- 
tures, we can estimate the visibility of the second or- 
der processes by comparing its contribution n 2 {co) to the 
contribution n-r(w). The result is 



n 2 (uj) 



1 



/ vj 2 {3\ k 2 [cosh(/3uj) — cos(ir/g)] 



n T (uj) ' ' (2?r) 1 /9+ 1 ^ a 



r(i-^-^) 



r(l-l/ 5 )sin(7r/ ff ) 

2 



(39) 



In order to eliminate the tunnelling amplitude 7 we can 
use the following relation with the total current / and 
the voltage 



ffr(l- l/g) sm{ir/g) 



ej u ( aeV 



1/9 



2n 2 ha \ hv 

Using this equation, we arrive at the following result 

n 2 (u>) 



(40) 



where 



MP,i,v) 



(2^) 1 ^ 1 /3fc 2 
2 ff r 2 (l-l/ ff )sin 2 (V.9)' 

h/3l ( hv \ 1/9 

aeV 



(41) 



(42) 



/iG3w) = [cosh(/3w) - cos(7r/ 5 )] (0w) 



(43) 
(44) 



r 1 



1 i/3u 
2g ~ ~2tT 



Additionally, we can determine the ratio of the second 
and fourth order tunnelling. Using l|24|) and <|29|) . one 

obtains 



n 2 (w) 2o , 
— ; r^r = — r(l - I/5) sm(7r/ 5 ) — 
714(0;) 7r K4 \ nv 



k 2 ( aeV\ ei 



7r < 45 > 




FIG. 3: Ratio of the secondary current to thermal broadening. 
In the inset: Ratio of the second order to the fourth order 
contribution 



In order to access information about the nature of cor- 
relations inside the host it is sufficient to compare the 
experimental data to the results in the lowest order in 
tunnelling, as the latter constitute the dominant contri- 
butions and are known in great detail. However, they are 
the more pronounced the closer to Ep j>{ui) is measured. 
In these regions the influence of both 774 and tit becomes 
important. Therefore, from the experimental point of 
view it is very important to possess detailed information 
about the energy window where only n 2 yields the ma- 
jor contribution. Having established the relations l|45|) 
and (|39|l between all three principal contributions to the 
high-energy tails we are now in a position to accomplish 
this task. 

The lower threshold of experimentally measurable cur- 
rents is expected to be about InA. Experimentally ap- 
plicable electric fields are in the range of lV//xm and 
the interface width is of the order of magnitude of lnm. 
Taking into account the amplification factor due to the 
sharpness of the emitter tip of ~ 10 3 we estimate the 
cut-off parameter V to be around IV. Assuming the 
experiment to be performed at T = IK, the resulting 
graphs for the ratios of n 2 {oj) to and 712(0;) to 

714(0;) are presented in Fig. [3] In order to compare the 
energy scales, we note that at a temperature of IK and 
a current of InA, fcgT r* 20IH/e. With our choice of pa- 
rameters, we arrive at the following values for the lower 
thresholds for g sw 0.5: n 2 > iit for u > 2ksT, and 
n 2 > 714 for u> > O.ObksT, so that at least at IK the 
processes of higher order in tunnelling can safely be ne- 
glected. The lower threshold for j>(ui) is then of the 
order of temperature. The upper threshold for the sec- 
ond order in tunnelling is given by to = eV as pointed 
out in Section ITll 
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V. CONCLUSIONS 

We have investigated field emission from strongly cor- 
related emitters with special focus on the high-energy 
tails of the energy resolved current. Using Keldysh dia- 
gram technique, we could isolate the relevant diagrams 
in fourth order in the tunnelling amplitude. Modelling 
the interaction in the host wire by a Luttingcr liquid 
with an interaction parameter g we could reproduce the 
already known leading asymptotics of the second order 
result, j> ~ uj 1 / 9 ^ 2 by an alternative approach. Using 
this method we succeeded in identifying the leading be- 
haviour of the high-energy tails in fourth order in the 
tunnelling amplitude, which turns out to be given by 
a; 1 ' 5-3 . Furthermore, we addressed the question of the 
observability of the secondary current by comparing the 



second order in tunnelling amplitude contribution with 
the broadening of the Fermi edge at finite temperatures 
as well as with the fourth order contribution. With typ- 
ical parameter values at hand we made an estimation 
of the energy window where only tunnelling processes of 
second order give dominant contributions to the particle 
energy distribution above the Fermi edge. 
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